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ABSTRACT 


Tangential mode, non-linear wave motion in a liquid propellant rocket 
engine is studied, using a two phase detonation wave as the reaction model 
Because the detonation wave is followed immediately by expansion waves, 
due to the side relief in the axial direction, it is a Chapman -Jouguet wave. 
The strength of this wave, which may be characterized by the pressure 
ratio across the wave, as well as the wave speed and the local wave 
Mach number, are related to design parameters such as the contraction 
ratio, chamber speed of sound, chamber diameter, propellant injection 
density and velocity, and the specific heat ratio of the burned gases. In 
addition, the distribution of flow properties along the injector face can 
be computed. Numerical calculations show favorable comparison with 
experimental findings. Finally, the effects of drop size are discussed 
and a simple criterion is found to set the lower limit of validity of this 
’’strong wave" analysis. 
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NOMENCLATURE 


area 

speed of sound 

function of Mach number, M 0 , defined by Eq. (5) 

specific heat at constant pressure 

diameter 

oxidizer to fuel mass flow ratio (= m 0 /m j) 
specific enthalpy 

total enthalpy in the wave fixed coordinate system 
thickness 

length and wave length 
mass flow rate 
Mach number 
pressure 
velocity 

heat release (= h^ - h°) 

mean radius of the annular chamber 

gas constant 

time 

temperature 

velocity component in the axial direction 
velocity component in the tangential direction 
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X. 
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y 

z 
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r 

p 

e 

v 

T 

Subscripts 

0 

1 

b 

c 

cr 

d 
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wave propagation velocity 
axial coordinate 
detonation length 
jet break-up length 
tangential coordinate (= R0) 
parameter defined in Eq. 11 
specific heat ratio 
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parameter = Vy 


y+1 

2(y-l) 


spacial density 
tangential coordinate 

non-dimensional tangential coordinate (= y/ x d ) 
period 

angular velocity 


condition in front of the wave 

condition behind the wave 

breakup 

chamber 

critical 

droplet 

fuel 
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o oxidizer 

£ liquid 
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s condition behind the leading shock 

Superscripts 

o condition at absolute zero 

bar ( — ) average value 
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I. INTRODUCTION 

Randomly distributed small explosions, sometimes described as 
M popping"\ often occur in liquid propellant rocket motors. Since the 
energy released per unit volume near the injector plate is high, these 
small explosions may often be amplified through detonation-like pro- 
cesses, and finally lead to steep fronted and self sustained pressure 

i 

2 

waves. This phenomenon is sometimes called "resonant combustion" . 
The pressure wave resembles the detonation wave in fundamental re- 
spects; it propagates supersonically in the circumferential direction at 
a nearly constant frequency, and the pressure ratio across the wave front 
has been recorded as high as thirty. 

The main objective of recent research in non-linear instability, 
theoretical as well as experimental (e.g. Refs. 2, 3, 4) is to under- 
stand the cause of the finite disturbances and their subsequent amplifica- 
tion. However, should the amplification not be successfully subdued, 
it becomes imperative that the engine designer know the effects of vari- 
ous design parameters on the final amplitude of the wave. Then the 
amplitude of the detonation -like wave can at least be reduced to a mini- 
mum strength in the design stage. Hence, in this analysis, it is assumed 
that the non-linear wave is fully developed at its final strength, a two 
phase detonation wave; the details of its evolution are set aside. 
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By studying the steady state problem, one can determine the weakening 

r 

and strengthening effects of various design parameters on the wave 
amplitude, and thus indicate the proper direction for a stable wave design. 
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H. PHYSICAL PICTURE OF ANALYTICAL MODEL 

There are several facts which have emerged from experimental 
studies of the strong wave instability studied here (e.g. Ref. 2); they 
lead directly to basic assumptions for the analytical model to be des- 
cribed. 

First, it is known that the waves rotate at a constant angular veloc- 
ity, and this is thus taken to be the case in the analysis. A somewhat 
subtle result of this observation is the condition of periodicity in the 
fluid properties. For example, since the wave velocity doesn't increase 
or decrease, the local properties immediately in front of and behind the 
wave are unchanged with time. In fact, at a given point on the chamber 
wall, an instrument would always measure the same distribution of pres 
sure (for example) with time, between and across each succeeding wave. 

Next, it is observed that the rotating waves are strong near the 
injector plate, but very weak near the nozzle, and that a wave doesn't 
’’wind up" on itself, that is that even though the strength of the wave 
varies in the axial direction, the wave angular velocity is the same at 
each axial location. This is not to say that the whole wave structure is 
in one plane, for it isn't; it simply means that the weak and strong parts 
of the wave each travel with the same angular velocity, no matter what 
the wave shape may be. This is also assumed in the model. 
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Finally, it is observed that the propellant mass flow rates are not 
significantly changed from their equilibrium running values when the ro- 
tating wave motion occurs, and this is also assumed in the model. 

The geometry of the combustion chamber under consideration is shown 
schematically in Fig. 1. In order to avoid the complexity arising from 
the radial variations in the governing equations, the analysis is limited 
to an annular chamber with a thin annulus. However, in a cylindrical 
chamber, because the wave travels tangentially at a constant angular 
velocity, it is only near the periphery of the chamber that the wave is 
strong enough to sustain reaction, and only near the injector plate that 
significant amounts of propellant droplets exist. Hence, there exists 
an annulus in which the assumptions made in the present analysis are 
clearly approximated. Just as in the familiar equilibrium rocket motor 
performance calculations, viscous forces are taken to be small com- 
pared to pressure forces and the heat transfered through the walls is 
considered to be small compared to the stagnation enthalpy of the gases 
in the chamber; thus, it is assumed that the chamber walls are frictionless 
and adiabatic. As a result, between waves, the gas flow is isentropic, 
and homoenergetic. 

A fully developed two phase detonation wave (or waves), rotating at 
a constant velocity, V (= HO), in the negative 6 direction, is considered. 

It is assumed that this wave front is planar and parallel to the axial direc- 
tion. It extends an axial distance x^ from the injector plate, where x^ 
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is the distance to which the propellants penetrate between successive 
waves. Since all propellants are consumed within x d , for x > x d> the 
wave is a shock rather than a detonation wave. In fact, as will be seen, 
the detonation branches into two oblique shocks at x = x d * For simplic- 
ity, the detonation is treated as a discontinuity, and all reaction is 
assumed to take place and be completed inside the wave. This assump- 
tion is valid as long as the time for droplet breakup and combustion 
behind the shock part of the detonation wave, is very small compared 
to the time between waves (or from another viewpoint as long as the 
wave thickness is small compared to the circumferential distance between 
waves). It is this condition which sets a limit on the analysis, as will 
be seen. 

Since the wave propagates supersonically, the period of the wave 

is extremely short (for example the typical wave period of an 11 in. 

5 6 

OD chamber is of the order of 470 nsec). It can be shown ’ for 
droplets of moderate size (for example 100 microns in diameter) which 
are typically found in present day rocket motor combustion chambers, 
that in one wave period the changes in droplet mass, momentum, and 
temperature are small compared to their initial values. Hence, it may 
be assumed that the droplets are mechanically and thermally frozen at 
their injected conditions in one wave period. In addition, the effects of 

5 

the droplets on the gas motion can be shown to be negligible . Hence, 
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to good approximation, the two phases near the injector plate may be 
considered as two independent systems, with phase change occurring 
only within the detonation wave. 

As a result of the side relief offered to the gases immediately behind 
the detonation wave, in the axial direction, the detonation wave is assumed 
to be a Chapman -Jouget wave. That is the velocity of the gases leaving 
the wave, relative to the wave, is sonic, so no more pressure pulses 
can travel upstream to weaken the wave further. 

Finally, it is assumed that the spray is dilute enough that droplet 
interactions may be ignored, that the gases follow the perfect gas law 
and are calorically perfect (constant specific heats), and that the energy 
available through reaction is not changed significantly by changing the 
mode of reaction. The first two of these assumptions are made for 
simplicity; the last is justifiable as long as the concentrations of the 
various species obtained at the end of wave reaction do not vary markedly 
from the corresponding values which the gaseous products would have 
in reaction taking place at the design condition. Certainly, as a first 
approximation, this variation in species concentrations may be neglected 
from the viewpoint that the propellant mixture ratio and hence the basic 
energy available is the same; only the mode of reaction is changed during 
transition from the design condition to the wave running condition. It 
should be noted that the main effects of different fuels and mixture ratios 
are introduced through this energy of reaction. 
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The analytical model described above has emerged as a result of 
efforts made to simplify the calculations, while still retaining and empha- 
sizing the essential mechanisms. At the nozzle end of the chamber, it 
resembles a pure acoustic model, because the waves are very weak waves 
which must travel therefore at essentially acoustic velocities for the 
mode considered, through gases which are essentially at conditions 
corresponding to those in the equilibrium running engine. This, plus 
the fact that all parts of the wave travel at the same angular velocity, 
explains the success of acoustic theory in predicting wave speeds. At 
the injector end of the chamber, however, the structure of the model 
has no relationship whatsoever to an acoustic model. Here, the waves, 
both detonations and shock waves are strong waves with pressure ratios 
as high as thirty. The properties of the flow field through which the 
waves travel are not even closely approximated by the properties in the 
equilibrium running engine. Typically, the high pressure and temperature 
immediately downstream of the detonation wave are rapidly decreased by 
strong expansions from values far above to values below the equilibrium 
running values. Hence the waves travel through "cold”, low pressure 
regions, compared to the design conditions; even though the velocity of 
the detonations is nearly the same as that of an acoustic wave in an undis- 
turbed combustion chamber, it is highly supersonic compared to the 
actual cooler gases through which it passes. Typical pressure distributions, 
illustrating this point, are shown later. 
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III. ANALYSIS 


The problem, in view of the rotating wave system, is basically unsteady 
in terms of a coordinate system fixed to the engine. This complexity can 
be removed by adopting a coordinate system fixed to the waves and thus rotat- 
ing at a constant tangential velocity, V . Such a system is illustrated in 
Fig. 2, where it is seen that in this case, the waves are stationary and the 
wall and droplets are moving at a velocity V . The periodicity found in the 
original fixed coordinate system is made manifest here by the repetitive 
nature of the flow picture; e. g. , a wave is seen each time a distance L 
has been traversed in the y direction in Fig. 2, with exactly the same 
conditions on either side of it as seen on the preceeding wave. Since no 
frictional forces or wall heat transfer is considered, the motion of the wall 
presents no computational problems. Because the two phases may be con- 
sidered as separate systems, the governing equations for the gases can 
be studied independently of the droplet motion. Finally, since radial vari- 
ations are neglected, the conservation equations for the gases may be 
finally reduced in this system to conventional two dimensional steady state 
equations. 

It is seen in Fig. 2, that the effects of the side relief are accounted 
for by an expansion fan which is centered at the triple point. The strength 
of the fan must be such as to bring the gas conditions back from their 
values immediately downstream of the detonation wave to the conditions 
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which existed immediately upstream of the detonation wave thus satisfy- 
ing the periodic condition. The strength of the oblique shocks must be 
such that the pressure must be continuous across, and the velocity vec- 
tors must be parallel on either side of, the slip line (the line which 
separates the flow which goes through the detonation wave from that which 
passes instead through the shock waves at the triple point). Downstream 
of the weak wave (characteristic) reflected at the intersection of the last 
wave of the expansion fan with the wall, the flow is parallel to the wall 
in the region bounded by the wall and the slip line : 

Governing Equations and Method of Solution 

As a result of the previously mentioned assumptions, it is easily 
5 

shown that the specific entropy, and total relative enthalpy which is 
defined as 

H = h + i(u 2 + v 2 -V w \ (1) 

must be constant along every streamline in the gaseous flow field between 
any two consecutive waves (see Fig. 2). However, since the detonation 
wave itself is one dimensional , there is no variation in the specific entropy 
and total relative enthalpy normal to the streamlines. Consequently, 
between the waves, the gaseous flow field is isentropic and homoenergetic. 
Furthermore, in order to satisfy the periodic condition, these two proper- 
ties must reach the same values in front of the wave in each and every 
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cycle Therefore, it is concluded that there should be no net change 
of the specific entropy and total relative enthalpy of the gases across 
the detonation wave itself. This in no way violates the Second Law of 
Thermodynamics, since the entropy of the mixture of gases and drop- 
lets does increase across the wave. 

The mixture entering the wave consists of burned gases and unburned 
propellant droplets, while the flow leaving the wave consists only of burned 
gases. The following non-dimensional quantities are defined: 



where the subscripts "o" and ”1” refer to conditions upstream and down- 
stream of the wave, respectively, and 0 denotes the energy released 
across the wave per unit mass of propellants; that is 

Q = h d - h° (3) 

5 

is the spacial density of the propellant droplets. It can be shown 
that there are eight equations (three from the jump conditions across 
the two phase detonation wave, two from the condition that the specific 
entropy and total relative enthalpy of the gaseous phase must be conserved 
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across the wave, two from the equations of state, and one from the 
Chapman -Jouguet condition) for the nine parameters associated with the 
detonation wave, i.e. P^, p 1 *, Tj*, V w , T q , P q , p Q , and v q . 
These equations may be written as explicit equations for the wave param- 
eters in terms of M q , the local wave Mach number, as follows: 


T 1 * = B 
p * = 




where B depends only on M q and y, 


(4a) 
(4b) 
(4 c) 
(4d) 

(4e) 

(4f) 

(4g) 

(4h) 
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B = 


(y + i) 


l 


1 + ^ n - 1 - M 


•3 


( 5 ) 


T q and p Q can be obtained from Eqs. (4f) and (4g) respectively, and v q 
and P q are calculated using the definition of M Q and the equation of state 
respectively; thus, 


v = M Vy R T~ 
o ooo 


P = p R T 
0 *0 o o 


(6a) 

(6b) 


It should be noted in Eq. (4f) that T c * is the dimensionless chamber 
temperature, T c /T q , and that the fact that Q* = T c * is the direct conse- 
quence of the fundamental assumption that the energy release is inde- 
pendent of the wave motion. Typical variations of T^*, P^*, V w *, T c *, 
p ( j*, and V w /a c with M q , are shown in Figs. 3, 4, and 5, for a representa 

tive value of y - 1. 26. In addition, P * = P /P and T /T , the pressure 

s so so 

and temperature ratios respectively across a shock wave at the same 
Mach number, M q , are shown for comparison with the corresponding 
ratios across the two phase detonation wave. 

Equations (4) do not afford a complete solution for the parameters 
of the problem, because, as noted previously, there are so far, only 
eight relations for nine unknowns. This is seen in Eqs. (4) by the fact 
that M q is needed to calculate the remaining parameters. The final 
necessary relation is found by considering the flow through the nozzle and 
relating this flow to that from the injector plate. 
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It is seen, in Fig. 2, that the flow field between the first slip line 

and the nozzle end of the chamber is very complicated due to the existence 

of oblique shock waves induced by the detonation wave. However, if these 

waves become sufficiently weak at the nozzle end, as assumed, the flow 

through the nozzle may be considered as a mean flow which satisfies the 

conventional one dimensional flow equations plus two dimensional small 

disturbances. To the degree of approximation used in this analysis, 

one may neglect these small disturbances, and employing a control 

volume which extends from the injector plate to the nozzle entrance and 

is bounded by the chamber walls, derive conservation equations (for 

mass, axial momentum, and total relative enthalpy) which relate conditions 

at the injector plate to the mean flow conditions at the nozzle end. Thus, 

the difficulties inherent in the irreversible flow region, where many 

oblique shocks with finite strength exist, can be avoided. As a result, 

and because the mass rate of propellant flow with waves is unchanged 

from its design, wave free value, one can relate the average pressure 

5 

on the injector plate to the design chamber pressure, P c ; that is, 



The average pressure on the injector plate may be written as 
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where t] (= y Q /x is the non-dimensional distance along the injector 
plate from the wave front to the point where the flow has expanded back 
to those conditions which exist upstream of the wave front, satisfying 
the periodic condition. The pressure distribution along the injector 
plate, P/Pj can be calculated using the method of characteristics 5 ’ ^ 
since the flow field between consecutive waves is irrotational. It is 
seen from Fig. 2 that the flow pattern resembles that of half of an under- 
expanded two-dimensional nozzle with the throat at the Chapman -Jouguet 
plane and centerline along the wall. The pressure distribution, P/P^, 
and Mach number, M, along the wall are plotted in Fig. 6. Finally, 

the pressure distribution can be integrated numerically to obtain the 

Vo 

average pressure; curves of r\ J (P/P 1 ) d? 7 , and J 0 are illustrated 

o 0 1 ^ 

in Figs. 3 and 4 for y = 1. 26. It should be noted that after P^ is obtained, 
a dimensional pressure distribution along the wall (pressure -space trace) 
can be obtained easily from Fig. 6. A typical wall pressure trace is 
shown later (Fig. 8). 
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From Eqs. (7) and (8) one can show that 



Also, since the total propellant mass flow, which comes from the injectors 


and passes through the (choked) nozzle, must also enter the detonation 
wave, 


V 


P„ X A 

c A d 

= p , V — 
d w c L 


o c 


Finally, Eqs. (9) and (10) may be combined to give, 


X, 

— J 
L J 1 


f 1 + r J 2 ) 


= z 


where and Z are defined as follows: 



( 10 ) 


( 11 ) 


(12a) 


(12b) 



Using Eqs. (4), one can write J, in terms of M q ; this function is shown 
in Fig. 4 for y - 1. 26. 

Equation (11) is the desired relation to give M q , as a function of y, 
Aj./A c , and x^/L. Since it is not possible to write an explicit functional 
form for M q , a numerical solution is necessary. A typical plot of a 
solution to Eq. (11) is shown in Fig. 7, where x^/L is plotted vs M q 
for various constant values of Z for y = 1. 26. Since y is a given value, 

Z varies only with the area ratio A^/A . The curves are found by setting 
Z equal to a given value, choosing a value for M q , calculating and 
(or using Fig. 4) and then calculating x d /L using Eq. (11). This is done 
for several values of M q , and repeated for several values of Z, leading 
to the curves shown in Fig. 7. 

Although a typical numerical calculation will be illustrated later, it 
is worthwhile to point out that at this stage of the calculations, all desired 
parameters and fluid property variations could be calculated if the follow- 
ing engine design parameter are known: 

(1) A/A c 

(2) y 

(3) P c 

(4) T c (a c ) 

(5) B 0 

(6) x d /L 
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where, again, P and T (or a ) are those values associated with the 
equilibrium running engine, and R q is the specific gas constant and thus 
equal to the universal gas constant divided by the molecular weight of 
the exhaust gases in the equilibrium running engine- The procedure is 
as follows: 

(1) Knowing A t /A and y, calculate Z, using Eq. (12b). 

(2) Knowing x^/L and Z find M q , using Fig. 7. This presupposes 
a plot such as Fig. 7 for the given y. 

(3) Knowing M q , calculate B (Eq.(5)) and then T^*, P^*, p ^*, v^*, 
V w *, T c *, P d *, and V w /a c , using Eqs. (4). 

(4) Knowing a and T , as well as V /a„ and T *, calculate V , 

' ■ c c w c c ’ w 

and T . 
o 

(5) Calculate P *, using Eq. (9), and knowing P , calculate P . 

C 0 0 

(6) Knowing P Q , T q , and R q , calculate p Q , using Eq. (6b). 

(7) Knowing P q , p Q , T q , and R q , as well as T^*, Pj*, p 1 *, Vj*, 
P d + , and M q calculate T^, P^, p^, v q , v^, and p^. 

The pressure distribution can be obtained by finding rj Q using Fig. 3, 
and then calculating y Q /L = V Q ’X^/ L. Then one can use Fig. 6, since P^ 
is known, to plot P vs y/L up to y Q /L, after which P = P q = constant up 
to y/L = 1 where a new wave exists. Thus, all quantities of interest can 
be obtained. It should be noted that knowing both R 0 and a c implies that 
one knows the fuel-oxidizer mixture ratio. 
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Calculation of Wave Width, 

In the previous section, it was shown that in order to obtain a general 
solution one must know not only the familiar design parameters of a 
rocket engine, but also x^/L. ft is this latter parameter which contains 
all the parameters associated with the injectors; one must be able to 
estimate this important parameter for various injector systems. 

is taken to be the furthermost penetration distance of spray for 
the given injection system in the time between successive waves; i.e. , 
it is assumed that no droplets exist beyond x^. For simplicity, it is 
assumed that the droplet distribution is uniform over x^, that is, that 
the wave is one dimensional. 

In the simple case where droplets are formed at the plane of the 
injector plate, x^ may be defined as the distance over which droplets 
move in one wave period, that is, 





a V 


c w 


(13) 


where u^ is the average droplet axial velocity, and would be considered to 
be known. However, it is noted that Eq. (13) generally underestimates 
the penetration of the spray for those injectors presently in use. Firstly, 
the spray tends to form at a finite distance from the injector plate, since 
a finite time is required for the liquid jets to be disintegrated. Secondly, 
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droplets in the spray are not uniform in size and thus are not uniformly 
distributed; larger drops penetrate further downstream. Hence, the 
penetration distance of the actual spray tends to be longer than the pene- 
tration distance calculated on the basis of the mean drop size. 

There is no known exact theory to predict the penetration distance 
of the spray during the rotating wave motion. The following analysis is 
an approximate calculation, wherein the difficulties of the non uniform 
distribution of droplets in the spray are ignored and a simple expression 
derived using the equations expressing the conservation of mass and axial 
momentum of the injected liquid jets. It is assumed that the liquid jets 
are disintegrated at some distance, x^, from the injector plate, and the 
jet breakup points are viewed as uniformly distributed sources, from 
which droplets issue. In the case where the jets are injected into the 
high speed transverse flow, x^ may be expressed by 


x i = u { ‘b 


(14) 


where u^ is the velocity of the liquid jet and t^ is the jet breakup time 

. 9 

given by 


d. 

t = 4 2 

b lq - q.l 


VT 


(15) 


In Eq. (15), d. is the jet diameter, q - q^ is the slip velocity between the 
jet and the surrounding gases, p^ is the density of the liquid jet, and p 
is the density of the surrounding gases. In the case of impinging jets, 
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the impingement distance may be used for x^, provided that the impinge- 
ment distance is less than the length calculated using Eq. (14). 

By using the conditions that the mass and axial momentum of droplets 
after breakup must be equal to the corresponding properties of the liquid 
jets coming into the point sources, a simple expression for x^/L can be 
derived^; 



u. 

V 


(16) 


w 

where u. is the mean jet velocity. In the case where bipropellants are 
used, then, 


u j = (TTf) u f + (TTf) u o (17) 

where f is the oxidizer to fuel ratio, and u f and u q are the axial velocity 

components of the fuel and oxidizer, respectively. Due to the difficulty in 

assessing the effects of the rotating wave motion on the jet breakup, x^ 

is approximated here by the predicted wave free jet breakup distance 

for the given injection system. 

Example Calculation; Comparison with Experiment 

A series of experiments on detonation- like wave motions were con- 
11-14 

ducted at JPL for various design conditions which are listed in 
Table I. In order to test the validity of the present analysis these design 
conditions were used to calculate various wave and flow properties for 
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comparison with experimental results. A typical calculation is presented 
here to illustrate the method of solution in detail. It should be noted that 
all the design conditions and experimental data presented in Table I are 
for cylindrical chambers; as mentioned previously, however, there is 
an annulus in the cylindrical chambers, where the assumptions used here 
are met. Only data from those transducers located in this annulus, i.e. 
near the corner made by the injector plate and the chamber wall, are 
used for comparison. Some are located on the injector plate and some on 
the wall. 

For the calculation to be made now, the design parameters of Ref. 11 

(Table 1) are used. In addition, y is chosen to be 1. 26. Thus, A./A = l/2, 

t c 

a = 3530 ft/sec, P = 300 psia, x. = 0. 75 in. , L = 34. 54 in. (for one wave), 

C C l 

f = 2. 8, q f = 148 ft/sec, q Q = 86 ft/sec, 9 f = 28° 18’ and 9 q = 15° 42’. Here 
q^ and q Q are the velocities of the fuel and oxidizer liquid injection jets 
respectively, and 9^ and 9 q are the corresponding angles made by these 
jets with respect to the axial direction. The injectors are of the imping- 
ing type and so xj is taken to be the axial impingement distance from the 
injector plate. 

From Eq. (12b), z is calculated to be 0. 2796. Next from Eq. (17), 

u. is calculated to be 93.6 ft/sec. Then, from Eq. (16), 

J 

x d /L = 0. 0432 + 0. 0265(a c /V w ) 
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(18) 



where V w /a c is given in terms of M q by Eq. (4h) . If Eq. (18) is plotted 

on Fig. 7 (dashed line) it is found to intersect the z = 0. 2796 curve at 

M = 2. 73, and x ,/L = 0. 0575. With this value of M , V /a =1. 847 

(from Fig. 3 or Eq. (4h)), P^* = 14. 74 (from Fig. 5 or Eq. (4b)), Jg = 18. 1 

(from Fig. 4, defined by Eq.(8b)), and P c * = 1.94 (Eq. (9)). Other desired 

parameters such as T^*, etc. , could be found in the same way by using 

Figs. 3 and 4 or Eqs. (4). Since P c = 300 psia, then P q = = 155 

psia, and P. = P /P i * = 2280 psia. Since a = 3530 ft/sec, V = a (V /a ) 
* [ 1 o 1 r c ’ w c w c 

= 6520 ft/sec. The calculated and experimental results are both shown in 
Table II for this and other cases. It is seen that the analytical results 
compare quite favorably with the experimental data, especially in the case 
of Ref. 11. 

After P^ is obtained, the dimensional pressure distribution along the 
injector plate can be obtained easily from Fig. 4, after determining ri 0 
from Fig. 3. In this case, for M = 2. 73, r\ = 4. 63, so y /L = 7) (x ,/L) 

= 0. 266. For y > y Q , of course, the pressure is constant at P . In 
Fig. 8, the calculated pressure distribution is shown and compared with 
an experimentally measured pressure trace obtained from Ref. 11. The 
agreement is seen to be excellent. 
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rv. EFFECT OF DESIGN CHANGES ON WAVE STRENGTH 

It is of interest to use the results given by the preceding analysis to 
predict the effects of various design parameters on the wave strength 
which is characterized here by the detonation wave pressure ratio, 

Pj* = P^/P . In each case, the effects are illustrated in the form of a 
plot of P^* versus the design parameter in question, other parameters 
being held constant at those values listed under Ref. 11 in Table I and 
used in the example problem, except where noted. 

1. Contraction Ratio, AjA„ 

v 

A plot of Pj* vs Aj./A c is shown in Fig. 9, for y = 1. 26 and for x^/L 

constant. A point may be calculated by choosing an M q , and finding the 

corresponding P^*, J^, Jg, and V Vf /^ c points from Figs. 3, 4, and 5. 

x^/L is here chosen to be the value calculated in the example calculations, 

x^/L = 0.0575. With x^/L, J^, and Jg known, z may be calculated using 

Eq. (11), and the corresponding A./A value may be found using Eq. (12b). 

t c 

By choosing several values of the parameter M q , the P^* and correspond- 
ing A./A values can be calculated and plotted as shown in Fig. 9. In addi- 
t c 

tion, a few points with x^/L varying as in Eq. (18) were calculated and are 
shown for comparison; this variation has a minimal effect. 

It is seen that as A^/A c increases, the pressure ratio, P^*, increases; 
the wave becomes stronger. Conversely, decreasing the contraction ratio 
has a stabilizing effect (wave becomes weaker). It can also be shown, 
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using Eqs. (4) that the wave speed increases or decreases the same as 
P^*. This result is also found experimentally^. 

2. Chamber Radius, R = L/2ir 

For a given x^, x^/L increases as R decreases, and it is seen in 
Fig. 10 that for other design parameters being held constant (most im- 
portantly a constant contraction ratio), decreasing R results in decreas- 
ing Pj*. Thus as engine size increases, for geometrically similar 
engines, the wave strength increases. Again, V w varies in the same 
manner as P *. 

The curve of P. * versus R can be calculated as follows. For 

A,/ A = 1/2 and y - 1. 26, z = 0. 2796. Then by choosing various values 
c c 

of M q , and finding the corresponding P^* and x d /L values from Figs. 5 
and 7, one can calculate L = x,/(x,/L) = 27rR for each M and thus for each 
Pi*. x d was chosen to be the value found in the example calculation, 
x d = 0. 0575 • 11 • it = 1. 99 in. (Note that Eq. (11) may be used rather than 
Fig. 7, with and given in Fig. 4.) 


3. Droplet Injection Velocity, u ^ 

In this case, the mass rate of flow of propellants is held constant 


and u. is varied by changing the injector areas. In addition, impinging 
1 

jets are assumed, so Eq. (16) is employed. It is seen in Fig. 11 that 
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as 11 ^ increases (which corresponds to an increase in x^), P^* decreases. 

The effect of increasing the jet velocity is to weaken the detonation wave. 

V varies in the same way as P„*. 
w J 1 

The computations for Fig. 11 are performed as follows. Equation 
(16) is rearranged to give, 



For these calculations, 2x./L = 0.0432 and a = 3530 ft/sec, the values 
from the example problem. Again, for A t /A c = 1/2 and y = 1. 26, z = 

0. 2796. The points are found by choosing values of M q , finding corres- 
ponding values of V w /a c , x^/L, and P^* from Figs. 3, 5, and 7, and 

calculating the corresponding value of u. from the above equation. 

J 


4. Chamber Pressure, P ., and Chamber Speed of Sound, a 

c ■ c 

P and a are considered simultaneously since they are intimately 
c c 

related by the mass flow relation, i. e. , m oc P /a_ - In the first case 
considered here m^ is held constant. In Fig. 12, P^* is plotted vs & c 
for the case of impinging jets, so Eq. (16) is used for x^/L; since m^ 
is constant, P varies also. P is plotted vs a in Fig. 13. This case 
corresponds to changing propellants or mixture ratios, holding other 


parameters constant. It is seen that as a c increases, P^* and P £ both 

increase. Moreover, since P *, plotted also in Fig. 12, remains 

c 
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relatively constant, then P q increases as P £ increases. Thus, not only 

the pressure ratio, P^*, but also the pressure levels (i.e. , P Q and P^) 

increase. Since V /a increases or decreases when P-,* increases or 

w c 1 

decreases, and P. * increases as a increases, V increases even more 

1 c w 

rapidly than a (i.e. , V = (V /a ) a ). 

* J c w V c c 

The curves in Figs. 12 and 13 may be computed as follows. From 
Eqs. (16), (9), (12b), and (10), one can show that 


u. a 
1 . c 

a = hi — 

c /x , 2x.\ V 


( x , ^x.\ 

T-t) 


w 


p * 

c 


■(’ 


J r Sy A c 


1 + L J 2/ Z r a. 


m A a 
p _c c 


c A c A t r ^ 


For these calculations, the values employed in the example calculations 
were used. Thus, 


A t /A c = l/2 y = 1. 26 


z = 0. 2796 


u. = 93.6 ft/sec 2x./L = 0.0432 m /A =1 lbm/in.^ sec 

] i p' c ' 


where the value of rhp/A c is that corresponding to the other parameters 
in the example calculation as seen in Column 1 of Table I. As in previous 
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calculations, various values of the parameter M q are chosen and corres- 
ponding values of x^/L, V w /a c , J and P^* are found from Figs. 7, 3, 4, 
and 5 respectively. Then the above equations may be used to calculate the 

corresponding values of a , P *, and P . 

c c c 

It should be noted that if an injection system where x^/L is a constant 

were considered, for this case where m^ is held constant, then from Fig. 7 

it is seen that since z is also a constant, M q would be a constant. Hence, 

from Fig. 5, P^ * would also have a constant value as a and P increased. 

However, since P c * (see Eqs. (9) and (12b)) would also remain constant, 

P would increase as P increased. Likewise P., =P *P * would in- 
o c 1 o 1 

crease as P c> The result would be a case where the pressure levels 

increased in such a way that the pressure ratio, P^*, remained constant. 

In addition since from Fig. 3, V w /a c would have a single value (for one 

value of M q ) , V w would vary directly as a £ . 

The second case to be considered here is that where a is held constant 

c 

and P c is varied, so that the propellant mass flow rate, nip, varies. How- 
ever, u. is held constant so that the m variation is due entirely to varia- 
J P 

tions in the cross sectional area of the injectors. Again impinging jets 
are assumed so that Eq. (16), holds, written again here as 


/x ■ 2x.\ 


2x.\ V 

w 

a 

a c c 
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Now, x.J L, u^, and a £ are constants, and it is seen from Figs. 7 and 3 

that since z is a constant, x,/L and V /a depend only on M . That is, 

the above equation could be written with M q as the only unknown parameter 

and thus it could be used to calculate the value of M for the case in ques- 

o 

tion. That is, M is a constant. Hence, from Fig. 5, P * is constant 
o 1 

and the plot of P^* vs P c would be simply a straight line, P^* = constant. 

Thus, from Figs. 3 and 7, or Eqs. (4), and from Eq. (9), V^/a^ x^/L, 

and P * have constant values. Then, as P increases;, with a constant, 
c c c 

P = P /P * increases and P„ = P * ■ P increases; the pressure levels 
o c c 1 1 o 

increase such that the pressure ratio remains unchanged. The wave speed, 
^w = <V a c> a c remains unchanged. 

5. Droplet Penetration Distance, 

The droplet penetration distance is a parameter controlled by the type of 

injector system employed and is thus indirectly a design parameter. In 

Fig. 10, the graph of P^* vs x^ indicates that as x^ increases, the wave 

strength decreases. In this case it can also be shown easily, from Figs. 7 

and 3 that as x , increases, V decreases. 

Calculations for the curve of P^* vs x^ in Fig. 10 are made as follows. 

A /A = 1/2 and y - 1. 26 are assumed so z = 0. 2796. Then, various values 
■ f c 

of M q are chosen and corresponding values of P^* and x^/L are found from 
Figs. 5 and 7. Finally x^ = (x^/L) L and L is chosen to be the value in the 
example calculation, L = 11 • n = 34. 54 in. 
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V. LIMIT OF THE ANALYSIS 


In Ref. 10, it is shown experimentally that an increase in drop size 
has a stabilizing effect, i.e. that the wave strength is decreased. On 
the basis of the present analysi s, this effect may be explained by noting 
that the larger droplets tend to penetrate farther downstream, leading 
to a larger x^, and hence a smaller P^*. However, the drop size effect 
is not formally incorporated into the present theory due to the funda- 
mental assumption that the wave is treated as a discontinuity. Now it 
15 

has been shown , that the reaction zone length in a two-phase detona- 
tion wave increases with drop size; hence losses due to lateral expansion 
occurring inside the reaction zone must become important at some point 
as the drop size is increased, and the thin wave assumption breaks down. 
Evidently the theory cannot be valid for droplets of all sizes. Although 
a detailed analysis of the structure of the reaction zone is not available, 
it is possible to set at least a rough limit on the drop sizes for which the 
analysis is valid. 

The structure of a spray detonation has been observed experiement- 
16 

ally . The study indicates that much of the mass transfer from the 
droplet is associated with a local explosion during droplet breakup. Hence, 
it is believed that a great portion of the energy released from the droplet 
may be associated with a blast wave initiating from the explosion site, 
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and that the interaction between the blast wave and the leading shock 
part of the detonation wave, may be mainly responsible for the energy 
transmission. On the basis of this observation, a simple criterion for 
the existence of a strong thin wave may be formulated. That is, the time 
required for a sound wave initiating at the explosion site of the disinte- 
grating drop to propagate upstream to the leading shock must be equal 
to or less than the time required for the same wave to propagate laterally 
the distance to which the detonation extends. If conditions immediately 

behind the shock are used to estimate the sound velocity and known ex- 

(17) 

pressions for droplet breakup times are used, this criterion finally 
leads to calculation of the critical drop diameter for a wave of given 
strength. A detailed derivation is given in Ref. (5); the resulting critical 
drop diameter is, 



where the subscript s refers to conditions behind the leading shock and 

m / p is the average volumetric flow rate, p being the average liquid 
P P P 

density is plotted in Fig. (14) for y - 1.26. If d is the average 
drop diameter calculated for the given injection system, then Eq. (19a) 
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leads to the following conclusions: 

(1) The strong wave analysis is valid if d < 'd . 

(2) Weak waves, which may not be properly described as detona- 
tions, are likely to occur if d > d cr> 

There are experimentally correlated equations to calculate d for dif- 

/1Q\ 

ferent types of atomization' For the case of impinging jets, for 
example, the following equation may be used to calculate the mass median 
drop diameter^) 

d = 6. 92 x 10 4 

where the diameters of the fuel and oxidizer injectors, d f and d Q 

respectively, are in inches and the jet velocities of the fuel and oxidizer, 

q^ and q Q respectively, are in ft/ sec. 

If d and d cr are calculated for each of the experimental conditions 

(5) 

covered in Table I, it can be shown' that only in the case of Ref. 11 is 
d less than d c ^; it is for this case also that the comparison between theory 
and experiment is most satisfactory. Evidently the waves in the remaining 
cases (see Table n) do not satisfy a strong wave criterion. 


, , 0.27 . 0.023] 
ld f d o 

0.74 0.33 

h % 


( 20 ) 
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VI. CONCLUSIONS 


This study is concerned with tangential mode detonation-like pheno- 
mena in liquid propellant rocket motors. It is found that the amplitude 
of the wave, measured by the pressure ratio across it (p^*), may be 

reduced by reducing the contraction ratio (A,/ A ), and chamber diameter, 

l c 

and by increasing the spray penetration distance during the wave running 
condition (x^). Expressions for x^ are provided. It is found that x^ 
increases as the jet velocity and projected jet breakup distance increase. 
In general, the variation of the wave speed follows the same trends as 
the wave pressure ratio. 

For a given contraction ratio, the effects of the chamber speed of 
sound and chamber pressure depend on the type of injector system used, 
and upon the propellant mass flow rate. In general, as a^ and p c increase 
the pressures immediately upstream and downstream of the wave, and the 
wave velocity, increase. In some cases, the pressure ratio across the 
wave also increase, and in others it remains unchanged. 

Finally, the limit of the strong wave analysis due to the finite drop 
size is discussed. A critical drop size, above which a strong wave 
probably cannot occur is derived, to set a lower limit of validity on the 
strong wave analysis. 
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0.173 

0. 101 
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B 956 
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Aver. 












B1093 


B1097 


Aver. 



V w 

(ftV sec) 

p i 

(psia) 

p o 

(psia) 

6116 

2392 

126 

6103 

2351 

121 

6116 

2469 

272 

6113 

2404 

176 

6520 

2280 

155 

6775 

882 

187 

6775 

1189 

185 

6775 

1036 

186 

7600 

2145 

151 

6871 

1145 

127 

6954 

998 

138 

6928 

1073 

133 

7700 

1830 

129 

6825 

877 

37 

6855 

891 

49 

6885 

455 

36 

6855 

741 

40.7 

7460 

1320 

44.6 







18. 98 473 


19.43 


9. 08 


13. 65 


14. 74 


4. 72 427 


6.43 427 


5. 57 427 


14.2 378 


9. 02 421 


7.23 416 


8.08 419 


14.2 374 


23.7 692 


18. 18 689 


12.64 686 






Table II. Comparison of experimental and analytical results. 
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Fig. 7. Xjj/ L versus M 0 for given values of z and y 
Dotted Lines Show Intersection of Eq. (18) with 
' Curve for z = 0.2796, giving M Q = 2. 73 for 
Example Calculation. 
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9. P^* versus A*/ A c for x^/ L = constant = 0. 0575 and 
T = 1. 26. The Points Indicated by •Show the Very 
Small Changes which Occur when an x^/ L which 
Varies According to Eq. (18) is employed. 
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R, X d (in) 


Fig. 10. Pj* versus and R for y = 1. 26 and A^/ A = 1/ 2 (z = 0. 2796). 

Calculations for were made Assuming a 5. 5 in. Radius 

Motor, as in Ref. 11, so that L = 34.54 in. Calculations 

for R were made Assuming x^ = 1.988 in. , corresponding 

to the numbers in the Example Calculation in the Test. 

(x,/L = 0.0575) 
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